In this paper, we present a well-organized strategy to estimate the FRACTIONAL advection-DIFFUSION EQUATIONs, which is an important class of EQUATIONs that arises in many application fields. Thus, Lagrange square interpolation is applied in the discretization of the FRACTIONAL temporal derivative, and the weighted and shifted Legendre polynomials as operators are exploited to discretize the spatial FRACTIONAL derivatives of the space-FRACTIONAL term in multi-termTIME FRACTIONAL advection-DIFFUSION model. The privilege of the numerical method is the orthogonality of Legendre polynomials and its operational matrices which reduces TIME computation and increases speed. A second-order implicit technique is given, and its stability and convergence are investigated. Finally, we propose three numerical examples to check the validity and numerical results to illustrate the precision and efficiency of the new approach.